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Abstract 

The surface detector array of the Pierre Auger Observatory is sensitive to Earth- 
skimming tau-neutrinos v T that interact in the Earth's crust. Tau leptons from v T 
charged-current interactions can emerge and decay in the atmosphere to produce a 
nearly horizontal shower with a significant electromagnetic component. The data 
collected between 1 January 2004 and 3 1 August 2007 are used to place an upper 
limit on the diffuse flux of v T at EeV energies. Assuming an E~' 2 differential 
energy spectrum the limit set at 90 % C.L. is El AN Vr /AE v < 1.3 x 10" 7 GeV 
cm -2 s _1 sr _1 in the energy range 2 x 10 17 eV < E v < 2 x 10 19 eV. 



The detection of Ultra High Energy (UHE) cosmic neutrinos at EeV (1 EeV = 
10 18 eV) energies and above is a long standing experimental challenge. Many ex- 
periments are searching for such neutrinos, and there are several ongoing efforts to 
construct dedicated experiments to detect them JHalzen e t al.(2002) Halzen(2007) 



|Falcke et al.(200 4) | . Their discovery would open a new window to the universe [Becker(2007) |, 
and provide an unique opportunity to test fundamental particle physics at energies 
well beyond current or planned accelerators. The observation of UHE Cosmic Rays 
(UHECRs) requires that there exist UHE cosmic neutrinos, even though the nature 
of the UHECR particles and their production mechanisms are still uncertain. All 
models of UHECR origin predict neutrino fluxes from the decay of charged pions 
which are produced either in interactions of the cosmic rays in their sources, or in 
their subsequent interactions with background radiation fields. For example, UHECR 
protons interacting with the Cosmic Microwave Background (CMB) give rise to the 
so-called 'cosmogenic' or GZK neutrinos | |Berezinsky et al.(1969)|. The recently re- 
ported suppression of the cosmic ray flux above ~4x 10 ly eV | Abbasi et al.(2007)| 
|Yam amoto(2 007)| Pierre Auger Collaboration(2007a)| as well as the observed corre- 
lation of the highest energy cosmic rays with relatively nearby extragalactic objects 
| Pierre Auger Collaboration(2007b)| both point to UHECR interactions on the infrared 
or microwave backgrounds during extragalactic propagation. These interactions must 
result in UHE neutrinos although their flux is somewhat uncertain since this depends on 
the primary UHECR composition and on the nature and cosmological evolution of the 
sources as well as on their spatial distribution [Engel et al.(2001) |Allard et al.(2 006)|. 

Tau neutrinos are suppressed in such production processes relative to v e or v^, 
because they are not an end product of the charged pion decay chain and far fewer 
are made through the production and decay of heavy flavours such as charm. Nev- 
ertheless, because of neutrino flavour mixing, the usual 1:2 ratio of v e to at pro- 
duction is altered to approximately equal fluxes for all flavours after travelling cos- 
mological distances [ Learned et al.(199~5)) . Soon after the discovery of neutrino os- 
cillations [Fuku da et al.(19 98)| it was shown that v T entering the Earth just below the 
horizon (Earth-skimming) [Fargion(2002)| |L~etessier-Selvon(2001)| [Feng et al.(2002)J 
can undergo charged-current interactions and produce r leptons. Since a r lepton can 
travel tens of kilometers in the Earth at EeV energies, it can emerge into the atmosphere 
and decay in flight producing an nearly horizontal extensive air shower (EAS) above 
the detector. In this way the effective target volume for neutrinos can be rather large. 

The Pierre Auger Observatory [Abraham et al.(2004)| has been designed to mea- 
sure UHECRs with unprecedented precision. Detection of UHECRs is being achieved 
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exploiting the two available techniques to detect EAS, namely, arrays of surface particle 
detectors and telescopes that detect fluorescence radiation. UHE particles such as pro- 
tons or heavier nuclei interact high in the atmosphere, producing showers that contain 
muons and an electromagnetic component of electrons, positrons and photons. This 
latter component reaches a maximum at an atmospheric depth of order 800 g cm~ 2 , 
after which it is gradually attenuated. Inclined showers that reach the ground after 
travelling through 2000 g cm" 2 or more of the atmosphere are dominated by muons 
arriving at the detector in a thin and flat shower front. 

The surface detector (SD) array of the Pierre Auger Observatory can be used to 
identify neutrino-induced showers [Capelle et al.(1998) Bertou et al.(2002), Zas(2005)|. 
The fluorescence detectors can also be used for neutrino searches | Aramo et al.(2005)l 
Miele et al.(2006)| but the nominal 10% duty cycle of the fluorescence technique re- 
duces the sensitivity. The electromagnetic component of neutrino-induced showers 
might reach the ground if the shower develops close enough to the detector, produc- 
ing a signal which has a longer time duration than for an inclined shower initiated 
by a nucleonic primary. Thus close examination of inclined showers enables showers 
developing near to the ground and those produced early in the atmosphere to be distin- 
guished. This allows the clean identification of showers induced by neutrinos, and in 
particular those induced by v T , with the SD JBilloir et al.(2007)| [Blanch Bigas(2007)| 
|Alvarez-Mufiiz(2007)l . 

Here we present the result of a search for deep, inclined, showers in the data col- 
lected with the SD of the Pierre Auger Observatory. Identification criteria have been 
developed to find EAS that are generated by r leptons emerging from the Earth. No 
candidates have been found in the data collected between 1 January 2004 and 3 1 Au- 
gust 2007 — equivalent to roughly one year of operation of the planned full array. 

The construction of the Southern Pierre Auger Observatory in Mendoza, Argentina, 
is currently close to being completed. It consists of an array of water Cherenkov tanks 
arranged in a hexagonal grid of 1 .5 km covering an area of 3000 km 2 that is overlooked 
by 24 fluorescence telescopes located at four sites around the perimeter. The array com- 
prises 1600 cylindrical tanks of 10 m 2 surface containing purified water, 1.2 m deep, 
each instrumented with 3 x 9" photomultiplier tubes sampled by 40 MHz Flash Analog 
Digital Converters (FADCs)[ Abrah am et al.(2 004)|. Each tank is regularly monitored 
and calibrated in units of Vertical Equivalent Muon ( VEM) corresponding to the signal 
produced by a fi traversing the tank vertically [Bertou et al.(2006)|. 

The procedure devised to identify neutrino candidate events within the data set is 
based on an end-to-end simulation of the whole process, from the interaction of the 
v T inside the Earth to the detection of the signals in the tanks. The first step is the 
calculation of the r flux emerging from the Earth. This is done using a simulation of 
the coupled interplay between the r and the v r fluxes through charged-current weak- 
interactions and r decay, taking into account also the energy losses due to neutral 
current interactions for both particles, and bremsstrahlung, pair production and nuclear 
interactions for the t lepton. The emerging r flux can be folded with the r decay 
probability to give the differential probability of r decaying in the atmosphere as a 
function of its energy and decay altitude, d 2 p T /dE T dh c . 

Modelling of the showers from r decays in the atmosphere is performed using 
the AIRES code |Sciutto(2002")) . The TAUOLA package QJadach etal.(1993)| is used 
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Figure 1 : FADC traces of stations at 1 km from the shower core for two real showers 
of 5 EeV. Top panel: electromagnetic component (9 ~ 22 °); bottom: muonic signal 
(6»-80°). 



to simulate r decay and obtain the secondary particles and their energies. Showers 
induced by the products of decaying ts with energies between 10 17 to 3 x 10 20 eV are 
simulated at zenith angles ranging between 90.1° and 95.9° and at an altitude of the 
decay point above the Pierre Auger Observatory in the range — 2500 m. Finally, to 
evaluate the response of the SD to such events, the particles reaching the ground in the 
simulation are stored and injected into a detailed simulation of the SD [Ghia(2007)|. 

A set of conditions has been designed and optimized to select showers induced by 
Earth-skimming v r , rejecting those induced by UHECR. The 25 ns time resolution of 
the FADC traces allows unambiguous distinction between the narrow signals induced 
by muons and the broad signals induced by the electromagnetic component (Figure[T|l. 
For this purpose we tag the tanks for which the main segment of the FADC trace has 
13 or more neighbouring bins over a threshold of 0.2 VEM, and for which the ratio of 
the integrated signal over the peak height exceeds 1 .4. A neutrino candidate is required 
to have over 60% of the triggered tanks satisfying these "young shower" conditions as 
well as fulfilling the central trigger condition [Abraham et al.(2004)| with these tanks. 
In addition the triggered tanks are required to have elongated patterns on the ground 



8 




length / width <speed> [m ns" 1 ] rms speed [m ns" 1 ] 



Figure 2: Distribution of discriminating variables for showers initiated by ts decay- 
ing in the atmosphere, generated by v T s with energies sampled from an E~ 2 flux 
(histogram), and for real events passing the "young shower" selection (points). Left: 
length/width ratio of the footprint of the shower on the ground; middle: average speed 
between pairs of stations; right: r.m.s. scatter of the speeds. See text for details. 

defining the azimuthal arrival direction (as expected for inclined events) by assigning a 
length and a width to the pattern and restricting its ratio (length/width>5). Finally, we 
calculate the apparent speed of the signal moving across the ground along the azimuthal 
direction, using the arrival times of the signals at ground and the projected distances 
between tanks. The average speed, as measured between pairs of triggered stations, is 
required to be compatible with that expected for an event traveling close to the hori- 
zontal direction by requiring it to be very close to the speed of light, in the range (0.29, 
0.31) m ns -1 with an r.m.s. scatter below 0.08 m ns -1 . These conditions are found 
to retain about 80 % of the simulated r showers triggering the SD. The final sample 
is expected to be free of background from UHECR-induced showers. In Figured we 
show the distributions of these discriminating variables for real events and simulated r 
showers. 

Over the period analyzed, no candidate events were found that fulfilled the selection 
criteria. Based on this, the Pierre Auger Observatory data can be used to place a limit 
on the diffuse flux of UHE v T . For this purpose the exposure of the detector must be 
evaluated. The total exposure is the time integral of the instantaneous aperture which 
has changed as the detector has grown while it was being constructed and set into 
operation. 

Calculation of the effective aperture for a fixed neutrino energy E v involves folding 
the aperture with the conversion probability and the identification efficiency. The iden- 
tification efficiency eft depends on the r energy E T , the altitude above ground of the 
central part of the shower h c (defined at 10 km after the decay point |Bertou et~a l.(2002) |), 
the position (x, y) of the shower in the surface S covered by the array, and the time t 
through the instantaneous configuration of the array. The expression for the exposure 
can be written as: 

r r Ev r°° d 2 » 

Exp = dn AE T dh c —^-B r , (1) 
Jn Jo Jq dE T dh c 

where f f 

' B T (E T ,h c ) = dt dxdy cos# eg [E T , h c , x, y, t] (2) 

where 9 and O are the zenith and solid angles. 
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The exposure is calculated using standard Monte Carlo techniques (MC) in two 
steps. The first integral deals with the detector-dependent part, including the time evo- 
lution of the array over the period T considered (eqf2]i. The integral in E T and h c 
involves only the differential conversion probability and B T (eqQ]i. The estimated sta- 
tistical uncertainty for the exposure is below 3%. 

The MC simulations require some physical quantities that have not been experi- 
mentally measured in the relevant energy range, namely the v interaction cross-section, 
the t energy loss, and the r polarisation. The main uncertainty in these comes from the 
QCD structure functions in the relevant kinematic range. We estimate the uncertainty in 
the exposure due to the v cross-section to be 15% based on the allowed range explored 
in | Anchordoqui et al.(2006)[ . The uncertainties in the t energy losses are dominated 



by the r photonuclear cross section. The 40% difference among existing calculations 
for the r energy losses JBugaev et al.(2004)| |Dutta et al.(2005)| |Aramo et al.(2005)| , 



which use different structure functions, is used as the systematic uncertainty. The two 
extreme cases of polarization give 30% difference in exposure and we take this as 
the corresponding uncertainty. The relevant range of the structure functions includes 
regions of Bjorken-x and squared 4-momentum transfer, Q 2 , where no experimental 
data exist. Only extrapolations that follow the behaviour observed in the regions with 
experimental data have been considered. 

We also take into account uncertainties coming from neglecting the topography 
around the site of the Pierre Auger Observatory | |Gora et al.(2007)| (18%). We adopt 
a 25% systematic uncertainty due to MC simulations of the EAS and the detector, 
dominated by differences between hadronic models (QGSJET | Kalmykov et al.(1997)) 
and SIBYLL lEngel et al.(1999)) ). 

Assuming a f(E v ) oc E~ 2 differential flux of v T we have obtained a 90% C.L. 
limit on the diffuse flux of UHE v T , whose level at 10 18 eV is representative for any 
smooth spectral shape: 

E 2 J(E V ) < l.dtojj x 10_7 GeV cm ~ 2 s_1 sr_1 ( 3 ) 



The central value is computed using the v cross-section from Ref. | Anchordoqui et al.(2006) | 



the parametrisation of the energy losses from Ref. [Dut ta et al.(2005)| and an uni- 
form random distribution for the r polarisation. The uncertainties correspond to the 
combinations of systematic uncertainties in the exposure as given above that lead to 
the highest/lowest neutrino event rate. The limit is applicable in the energy range 
2 x 10 17 — 2 x 10 19 eV, with a systematic uncertainty of about 15%, over which 
90% of the events are expected for f{E v ) oc E~ 2 . In Figure [3] we show our limit 
adopting the most pessimistic scenario for systematic uncertainties. It improves by a 
factor ~ 3 for the most optimistic one. For energies above 10 20 eV, limits are usually 
quoted as 2.3 /Exp x E v for different energy values (differential format), while at lower 
energies they are usually given assuming an E~ 2 flux (integrated format). We plot the 
differential format to demonstrate explicitly that the sensitivity of the Pierre Auger Ob- 
servatory to Earth-skimming v r peaks in a narrow energy range close to where the 
GZK neutrinos are expected. 

The Earth-skimming technique used with data collected at the surface detector ar- 
ray of the Southern Pierre Auger Observatory, provide at present the most sensitive 
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Neutrino Energy [eV] 



Figure 3: Limits at 90% C.L. for a diffuse flux of v T from the Pierre Auger Observa 
tory. Limits from othe r experiments |Ach terberg(2007) |Ackermann et al.(2007) 



Martens(2007) 



Aynutdinov et al.(2006) 



Barwick et al.(2006)| |Gorham et al.(2004)| |Lehtinen et al.(2004)) 



Kravchenko et al.(2006) 



are converted 

to a single flavour assuming a 1 : 1 : 1 ratio of the 3 neutrino flavours and scaled 
to 90% C.L. where needed. Two different formats are used: differential (squares) 
and integrated (constant lines). The shaded curve shows the range of expected 
fluxes of GZK neutrinos from Ref. JEngel etal.(2001)| [Allard et al.(2006)| , although 
predictions almost 1 order of magnitude lower and higher exist. 

bound on neutrinos at EeV energies. This is the most relevant energy to explore the 
predicted fluxes of GZK neutrinos. The Pierre Auger Observatory will continue to take 
data for about 20 years over which time the limit should improve by over an order of 
magnitude if no neutrino candidate is found. 
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